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We provide further characterization of non-Markovian quantum dynamics based on the concept 
of divisible dynamical maps. In analogy to entanglement witness we propose a non-Markovianity 
witness and introduce the corresponding measure of non-Markovianity. We also provide characteri- 
zation of non-Markovianity in terms of entropic quantities, fidelity and Wigner-Yanase-Dyson skew 
information. 
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I. INTRODUCTION 



if Lt has the standard form [H, Q for all t > 0, that is, 



The dynamics of open quantum systems attracts nowa- 
days considerable attention. It is relevant not only for a 
better understanding of quantum theory but it is funda- 
mental in various modern applications of quantum me- 
chanics. Since the system environment interaction causes 
I— ]■ dissipation, decay and decoherence it is clear that the dy- 
] namics of open systems is fundamental in modern quan- 
tum technologies, such as quantum communication, cryp- 
tography, computation and quantum metrology. 

The traditional approach to the dynamics of an open 
quantum system consists in applying a suitable Born- 
Markov approximation leading to the celebrated quan- 
tum Markov semigroup [H, Q which neglects all mem- 
ory effects. Recent theoretical activity and technological 
progress show the importance of more refine approach 
based on non-Markovian evolution. Non-Markovian 
quantum dynamics becomes in recent years very active 
field of both theoretical and experimental research and 
there are a lot of papers devoted to this topic ( see e.g. 
Q-fl^ and references therein). 

Surprisingly, the concept of (non)Markovianity is not 
uniquely defined. One approach is based on the idea of 
the composition law which is essentially equivalent to the 
idea of divisibility [l^ . A dynamical map At is divisible if 
At = Vt.sAs and Vt^s is completely positive and trace pre- 
serving for all t > s, that is, it gives rise to 2-parameter 
family of legitimate propagators. The essential property 
of Vf.s is the following (inhomogeneus) composition law 



Vt,u 



(1) 



for all < > s > M. It is clear that ([T]) generalizes semigroup 
property. This app roach was used by Rivas, Huelga and 
Plenio (RHP) [20| to construct the corresponding mea- 
sure of non-Markovianity which measures the deviation 
from divisibility. In this paper we assume that Marko- 
vian dynamics is represented by divisible dynamical map. 
It should be stressed that Markovian dynamics (divisible 
map) is entirely characterized by the properties of the lo- 
cal in time generators Lt, that is, if At satisfies At — LfAt, 
then At corresponds to Markovian dynamics if and only 



with time dependent Hamiltonian Hf and noise operators 
14 (t). 

A different approach is advocated by Brcuer, Laine 
and Piilo (BLP) in Ref. [iH. BLP define non-Markovian 
dynamics as a time evolution for the open system char- 
acterized by a temporary flow of information from the 
environment back into the system and manifests itself as 
an increase in the distinguishability of pairs of evolving 
quantum states: 



o-(pi,P2;t) 



1 d 
2di 



|At(pi - P2)||l 



(2) 



where = Tr\/ A^A denotes the trace norm. Accord- 

ing to [l9| the dynamics At is markovian iff cr(pi, P2'^t) ^ 
for all pairs of states pi, p2 and t > 0. Optimizing over 
pi , p2 enables to construct suitable non-Markovianity 
measure [l^ (see the recent paper [2l| discussing the 
properties of the optimal pair pi,p2)- It turns out that 
condition <7{pi, p2;t) < is less restrictive than require- 
ment of complete positivity for Vt^s and one can con- 
struct At which is non-Markovian (not divisible) but still 
gives rise to the negative flow of information (see (22| - [23 | ) . 
Other measures of non-Markovianity based on quantum 
Fisher information (sOj, fidelity [3l|, departure from di- 
visibility 27 1 and quantum mutual information [28| , were 
proposed as well. 

In the present paper we provide further characteriza- 
tion of non-Markovian quantum dynamics. We stress in 
our approach Markovianity of evolution corresponds to 
the divisibility of the corresponding dynamical map. In 
the next section we propose a non-Markovianity witness. 
In analogy to entanglement witness non-Markovianity 
witness is defined by Hermitian not positive operator 
in T-L^T-L. Using this concept we introduce a measure 
of non-Markovianity which is essentially equivalent to 
the RHP measure. We also provide characterization of 
non-Markovianity in terms of entropic quantities, fidelity 
and Wigner-Yanase-Dyson skew information. Simple ex- 
amples illustrate the differences and similarities between 
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these characteristics. In particular we provide another 
example of quantum evolution which is non-Markovian 
(not divisible) but still satisfies BLP criterion ([2]). Final 
conclusions are collected in the last section. 



II. NON-MARKOVIANITY WITNESS 

Let us recall, that if 8 : T('H) T{T-L) is a linear trace 
preserving map, then £ is a quantum channel if and only 
if 



(]l®£)X||i < ||X||i 



(3) 



for all X = e Bin^n). Traditionally TCH) denotes 
the vector space of trace class operators, i.e. x G T('H) 
if = TrV xx^ < oo. It is clear that if H is finite 

dimensional, then T{H) = B{T-L). Actually, is a 

contraction for all X not necessarily hermitian. Hence, a 
dynamical map At is Markovian if and only if 



\t{X) := -\\{t® kt)X\\i < 



(4) 



Definition 1 We call X^ ^ X non-Markovianity wit- 
ness for At iff At (A) > for some t > 0. 

Note, that if A > then {t® kt)X > and hence 
11(1® At)A||i = TrA which implies At(A) = 0. There- 
fore, a necessary condition for A to be non-Markovianity 
witness is A ^ 0. Recall, that it is also a necessary 
condition for an entanglement witness. 

We can propose a natural measure of non- 
Markovianity 



AA[At]- sup / \t{X)dt 

\\X\\i = l J\t{X)>0 



(5) 



that is, the formula ([5]) choses the optimal witness and 
measures the violation of At (A) < along the trajectory. 
Note, that condition Q implies 



^l|Ata:||i <0, (6) 
for all X = E 13(71), and ([6|) implies BLP condition 



^||At(p-a)||i<0, (7) 

for all density operators p and a in %. It is, therefore 
clear, that BLP definition of non-Markovianity is less re- 
strictive. 

Example 1 Consider pure decoherence of a qubit sys- 
tem described by the following local generator 



LtP = ]p{t){crzP(Tz - p) 



(8) 



The corresponding evolution of the density matrix reads 

-r(t) 



Pt 



Pii Pi2e 
pi2e"'"(*) P22 



(9) 



where T{t) ~ j{T)dT. The evolution is Markovian 
iff 7t > 0. Taking Aq = ^(Jx®<^x one finds Af(Ao) = 
-7(t)e"^^*). Hence Af(Ao) > whenever -f{t) < 0. We 
claim that Aq optimizes ([5]). Let us recall, that RHP 
compute 

m . hm + , (10) 

e 0+ e 

and get g{t) = \ j{t)\ whenever j{t) < 0. One has 

A/RHpiAt] ^^AFfc , (11) 

k 

where AFfe = |F(ifc -I- Ak) - F(tfe)| and j{t) < for t e 
{tk,tk + Afc). Similarly, one finds 



-r(tfc+Afc) _ p-r(tfc)| 



(12) 



It is clear that 7V[At] > if and only if A/RHpiAf] > 0. 

Suppose that a dynamical map At possess time- 
independent eigenvector, that is. At/ = p,tf, where /xt 
belongs to the unit disc in the complex plane. Markovian- 
ity implies that ^|/J.t| < for alH > 0. There are many 
examples of quantum dynamics where the off-diagonal 
elements behave according to 



M\^){:j\) = G,,m){J\ , 



(13) 



Hence, if ^|Gij(t)| ^ for at least one pair (ij), then Af 
is non-Markovian. 

Note, that if all eigenvectors fa of At are time inde- 
pendent 



^tfa = Pa{t)fa , Q! = 1 , . . . , (dimH)^ 



(14) 



then AtA„ = A„At, for all t,u > 0. One may call it 
commutative dynamics. If the dynamics is Markovian, 
then 

|l|At/o||i = |K(t)|||/a||i<0. (15) 

Therefore, for a class of commutative dynamical maps 
Markovianity implies monotonicity of all |/iQ,(i)|. Actu- 
ally, pure decoherence dynamics belongs to this class. 
One has 

M\k){k\) = \k){k\ (k = 1,2) , At(|l)(2|) = e-^W\l){2\ , 

and hence /ii(t) — p-2{t) = 1 and ^^{t) — iii{t) = e^'"^*-'. 
Hence Markovianity implies ^{t) > 0. 
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Example 2 Consider the dynamics governed by the lo- 
cal in time generator 



LtP = (ujtTrp- p) , 



(16) 



where ujt is a family of Hermitian operators satisfying 
Trcjt = 1. The above generator gives rise to Markovian 
evolution iff Lf has the standard form 0,0 for alH > 0, 
that is, iff -f{t) > and wt defines a legitimate state, 
i.e. Wt > 0. The corresponding solution of the Master 
equation pt = Ltpt with an initial condition pt=Q = p 
reads as follows 



Pt 



V+[l-e-^(*)]f)tTrp, 



(17) 



where fit = ^r(t}_i Jq li''')^^^'^'' ^rdr (note, that if p is 
density matrix then Trp = 1). It is therefore clear that Lt 
generates a legitimate quantum evolution iff T{t) > and 
^{t) > 0, that is, fif defines a legitimate state (note, that 
Trfif = 1). In particular, if = a; is time independent, 
then 17t = w and the solution simplifies to a convex com- 
bination of the initial state p and the asymptotic invari- 
antstatew: = e~^(*V+[l — e^^(*)] w. One easily shows 
that the evolution is Markovian iff 7(t) > and uJt is a le- 
gitimate density operator (that is, Wt > 0). Consider now 
the BLP condition One has pt ^ (Jt — e~'"(*'(p — a) 
and hence 



d 
di 



\Pt - o-tlli 



-7(t)e- 



|p- cr||i < 



implies only j{t) > but says nothing about uJt- It 
shows that any which gives rise to fit > leads to 
the evolution satisfying condition ^ but only > 
gives rise to Markovian dynamics. Hence, we may have 
non-Markovian dynamics (wt ^ but fit > 0) which 
satisfies BLP condition ^ for all t > 0. Clearly, 
such non-Markovian dynamics has vanishing BLP non- 
Markovianity measure. 



III. ENTROPIC WITNESSES 



for each p. In particular, if ctq = I/n is maximally mixed, 
then S{p II (To) = '^ogn — S{p), and hence the formula ([2T|) 
reduces to 



dt 



SiAtip)) > , 



(22) 



that is, the von Neumann entropy monotonically in- 
creases for each initial state p. 

Consider once more pure decoherence of a qubit form 
Example [T] Note, that Ltl ~ and hence the maximally 
mixed state is invariant. Therefore, Markovianity implies 
(22). Now, for the 2-level system 



(23) 



where > are eigenvalues of pt- Hence ^ S{pt) > 
ii < 0. One easily finds 



- l±V(pn-p22)2 + |pi2Pe-2r(t; 



It is therefore clear that S{pt) monotonically increases if 
and only ii t{t) = ■y{t) > 0. 

The above scheme may be immediately repeated for 
the well known families of generalized Renyi Sa and Tsal- 
lis Tq relative entropies 



SaipW C) 



1 



a — 1 

for a e [0, 1) U (l,oo), and 
Tqip\\cr) = 



log 



Trp^o-i-" 



l-TrpV^-? 



1-g 

for g G [0, 1). Note, that in the limit 



(24) 



(25) 



Let us recall that the relative entropy defined by 

S{p\\ a) = Tr{p[\ogp^ log a]) . (18) 

(one assumes that S{p \ \a) = oo when supports of p and 
a do not satisfy suppp C supper) enjoys 



S{£{p)\\£{a))<S{p\\a) 



(19) 



for any quantum channel £. Therefore, if At is a dynam- 
ical map, then S{At{p) \\ Af((T)) < S{p \\ cr) for any t > 0. 
Hence, if At corresponds too Markovian evolution, then 



dt 



5(At(p)||At(a)) <0 



(20) 



for each pair of initial states p and cr. Note, that if (Tq is 
an invariant state, i.e. At(cro) = ctq, then (|20l) simplifies 
to 



dt 



^(At(p)lko) <0 



lim (p 1 1 cr) = lim T„ (p 1 1 cr) = S'(p 1 1 cr) , 

CK — > 1 Q — > 1 

one recovers von Neumann relative entropy. It turns out 
d, that if is a quantum channel then Sa satisfies 
S'a(£:(p) ||£'(cr)) < Sa{p\\cr) foi-ae [0,1)U(1,2] and the 
same applies for Tq for g e [0, 1). If At is a divisible map, 
then 

5„(Atp||At(T) <0 , ^T,(Atp||Ata) <0 , (26) 



dt 



for a e [0, 1) U (1, 2] and q S [0, 1). Again, if the maxi- 
mally mixed state is invariant, then (j26p gives rise to 



dt 



SaiMp)) > 



d 

di 



TqiAtip)) > , (27) 



(21) which generalize (P^ . 
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IV. FIDELITY WITNESS 

Given two density operators p and a one defines 
Uhlmann fidelity 

2 



F{p,<7) = Tr JsfpcJ^ 



(28) 



Equivalently, one lias F{p,a) — \\^/p^/^\\l which shows 
that F{p, tr) = F((t, p). One proves 



1 - F{p, a) < D[p, a] < ^l-F{p,aY 



(29) 



Moreover, if £ is a quantum channel, then 
F {£ (p) , £ {(t)) > F{p,a) which implies that for the 
Markovian evolution one has 

|f(A,(p),A,(ct))>0. (30) 

Again, if ctq is an invariant state, then Markovianity im- 
plies ^ 1 1 -y/iot V^l 1 1 — 0- In particular, if cto is maximally 
mixed, then 



^IIVP^Ik = |Try^>0. 

Note, that the above condition is equivalent to 
d 



dt 



(31) 



(32) 



If tTo ~ |V'o)('0o| is pure, then ([50)) reduces to a very 
simple condition 



(33) 



for all t > 0. It means that the overlap of pt with an 
invariant vector state j-^o) monotonically increases. 

Example 3 (Spin-boson model) Consider once more 
the evolution of a qubit system described by the following 
local generator 



Ltp 



is{t) 



1 



[cr+cr_,p] +7(t)((T_pcr+ - -{cr+a^,p}) 



where s(t) = — 2Im and j{t) - 
hmction G{t) satisfies non-local equation 



(34) 

2Re , and the 



G{t)^- / f{t-T)G{T)dT , G(0) = 1 , 

Jo 

with f{t) being a correlation function of the (bosonic) 
reservoir. The standard raising and lowering operators 
read: a+ = |2)(1| and (t_ = |1)(2|. The corresponding 
evolution is given by the following formulae 

Pll(i) = Pll + (1 - \G{t)\^)p22 , P22{t) = \G{t)\^p22 , 

and the off-diagonal element pi2{t) = G*{t)pi2- It is 
clear that the ground state \^q) = |1) defines an invariant 
state, and hence Markovianity implies 

^(V'o|pt|V'o) = ^Pii(t)>0, (35) 

which is equivalent to ^ |G(i)p < and hence to "f[t) > 
0. 



V. WIGNER-YANASE-DYSON SKEW 
INFORMATION 

Consider the following quantity 

i{p,x)^-]p:v[^p,xY , 



(36) 



where X'^ = X is an observable, introduced by Wigner 
and Yanase [s^. One finds I{p,X) = Tr[pX2 _ 
■y/pATy/pX] and hence if p is pure then /(p, X) reduces to 
the variance V{p,X) Tr(pX2) - Tr(pX)2. This quan- 
tity was generalized by Dyson 

/p(p,^) = -^Tr[p^X][pl-^X] , (37) 

for arbitrary < p < 1, and its is called Wigner- Yanase- 
Dyson skew information. The convexity of / = I1/2 was 
already proved by Wigner and Yanase [H], and for the 
general case p S (0, 1) ~ the celebrated Wigner- Yanase- 
Dyson conjecture - it was proved by Lieb [3; 

It was proved by Petz 
channel, then 



s proved 



that if f is a quantum 



lp{£{p),X)>Ip{p,£*{X)), (38) 

where £* denotes a dual channel (Heisenberg picture) 
defined by Tv{£{p)X) = Tr(p£*(A:)). Hence, if ctq is an 
invariant state of A(, then 

/p(ao,X)>/p(ao,A*(X)) . (39) 

Similarly, if Xq is an invariant observable (constant of 
motion), i.e. Aj (A"o) = Xq, then 

/p(At(p),Xo)>/p(p,Xo) . 

In conclusion, if At is Markovian, then 



(40) 



-/p(ao,A:(X))<0 



and 



dt 



Ip{At{p),Xo) > . 



(41) 



(42) 



(V'ol^t'l^o)-(V'ol^tl^o) 



< 



(43) 



Interestingly, if cto = |V'o)(V'ol is a pure state, then Marko- 
vianity implies 

d_ 

di 

where Xt — A^ (AT). It shows that for Markovian dynam- 
ics dispersion in the invariant pure state monotonically 
decreases. 

Consider spin-boson model defined in Example 2. One 
finds for the evolution in the Heisenberg picture 

Xu(t) = Xu , X22{t) = (1 - |G(i)P)Xn + |G(t)|2X22 , 

and the off-diagonal element Xi2{t) = G{t)Xi2. Taking 
as an invariant state (Tq = the formula (j43p implies 

d 



dt 



[Xn + \Git)\^\Xi2\^] <0 , 



(44) 



which is equivalent to ^ |G(t)p < and hence to 7(t) > 
0. 
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VI. CONCLUSIONS 

We have provided several criteria enabling one to wit- 
ness the non-Markovianity of quantum evolution. Note, 
that passing to the Hcisenberg picture (that is, using the 
dual map A^) one may reformulate the formula Q as 
follows: A( corresponds to the Markovian dynamics in 
the Hcisenberg picture iff 

A,(X):^|||(]1®A*)X||<0, (45) 

where ||^|| denotes the operator norm in B{H(E)'H). 
Both criteria are perfectly equivalent: Xt{A) < for all 
A € T{n (g) and t > if and only if At (B) < for all 
B e B{H(^H) and t>0. Interestingly, Wigner-Yanase- 
Dyson skew information merges both pictures and pro- 
vides the constraints for Markovianity of quantum evolu- 
tion given by (pij) and (|42p . It is clear that the witnesses 
of non-Markovianity presented in this paper may be im- 
mediately generalized. 

One may for example try to use generalized entropic 



measures [1^ Sf{p) = Tr/(p), where f{p) is a smooth 
strictly concave real function defined for p G [0, 1] satisfy- 
ing /(O) — /(I) = 0. They provide further generalization 
of Rcnyi and Tsallis entropies. 

Note, that if _F is a function satisfying monotonicity 
condition F{£[p)) < F{p) for any quantum channel 5, 
then Markovianity of At implies ■^F{At{p)) < for all 

t ^0. Similarly, if F{[lA(g> SKpab)) < F{pab), then 
^F([llA«iAt](pAB)) < for Markovian evolution. This 
property was used by Rivas et. al. by taking as F 
genuine entanglement measure. Similarly Luo et. al. (28j 
take as F the mutual information. Fidelity and relative 
entropy may be replaced by any function G satisfying 
monotonicity condition G'(£'(p), £(cr)) <G{p,(t). Marko- 
vianity of At implies j-fi{Kt{p), ht{.o-)) < for aU t > 0. 
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